Tunneling broadening of vibrational sidebands in molecular transistors 



m 
O 
o 

(N 
Or 



Karsten Flensberg 
Laboratory of atomic and solid state physics, Cornell University, Ithaca, 
NY 14853, USA and 0rsted Laboratory, Niels Bohr Institute fAPG, 
Universitetsparken 5, 2100 Copenhagen, Denmark. 
(Dated: February 6, 2008) 

Transport through molecular quantum dots coupled to a single vibration mode is studied in the 
case with strong coupling to the leads. We use an expansion in the correlation between electrons 
on the molecule and electrons in the leads and show that the tunneling broadening is strongly 
suppressed by the combination of the Pauli principle and the quantization of the oscillator. As a 
consequence the first Frank-Condon step is sharper than the higher order ones, and its width, when 
compared to the bare tunneling strength, is reduced by the overlap between the groundstates of the 
displaced and the non-displaced oscillator. 
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I. INTRODUCTION 

In recent years, the field of molecular electronics has 
renewed the interest in transport through mesoscopic sys- 
tems with strong electron-phonon coupling. A number of 
experiments have been reported to demonstrate trans- 
port through single molecules some of which 
show signs of vibrational sidebands. In the pioneering 
experiments by Park et ali it was shown that electron 
transport through a single 6*60 molecule was strongly in- 
fluenced by a single vibrational mode. The single phonon 
mode was associated with the motion of the molecule 
in the confining potential created by the van der Waals 
interactions with the electrodes. Later similar devices 
with more complicate molecules were investigated^ and 
they also showed excitation spectra which possibly could 
be associated with emission of vibrational quanta, so- 
called Frank Condon peaks in the differential conduc- 
tance. Some of these devices^ furthermore exhibited a 
peak stemming from the Kondo effect proving that the 
tunnel coupling to the leads was rather strong. 

Theoretically there has been a large amount of work 
on the problem of tunneling through a single level with 
coupling to phonon modes. The approaches fall into two 
categories. The first category is the kinetic equation ap- 
proach, which is relevant in the weak tunneling limit. 7,8 
In the kinetic equation approach, it is essential that the 
excited vibrational levels are allowed to relax through 
coupling to a bath. For a large phonon relaxation rate 
one can assume an equilibrium phonon distribution, oth- 
erwise a non-equilibrium distribution function must be 
determined from the kinetic equations^ The coupling to 
the dissipative environment of the molecule leads to ad- 
ditional broadening of phonon sidebands. 8 

The second category, in which also this paper falls, 
deals with the strong tunneling limit. The first ap- 
proaches were an exact solution of a simplified situation, 
where only one electron is present)2ilfl i.e., the presence 
of the Fermi sea is ignored. The result of this exact solu- 
tion is that the phonon sidebands in the tunneling density 
of states all are Lorentzians with the same width. The 
exact solution amounts to decoupling the electron and 
oscillator displacement operators. This decoupling ap- 



proximation has been used in a number of papers^iS^ 
Below we refer to this approximation as the "single parti- 
cle approximation" (SPA) and we demonstrate explicitly 
that this approximation is only valid for high energies, 
i.e., for electrons or holes sufficiently far from the Fermi 
surface, or in the weak tunneling limit where the rate 
equation is valid. For the many-body system many de- 
cay channels are in fact blocked by the Pauli principle. 
Other many-body works that include Fermi sea effects 
have been done^i^*^ and the blocking of the line width 
has also been alluded toi 14 i 15 

In this paper, we focus on the broadening of the Frank- 
Condon peaks in the differential conductance due to tun- 
neling broadening, but in the non-Kondo regime. While 
this case cannot be solved exactly, an approximation, 
which is believed to be valid above the Kondo tempera- 
ture and which includes the Pauli principle, but ignores 
correlations between lead electrons, is developed. The 
main result is that the combination of the Pauli principle 
and the splitting of the spectral weight into phonon side- 
bands severely limits the phase space for tunnel broad- 
ening. Hence the peaks are much narrower than the bare 
tunneling rate would suggest. We find a simple analytic 
result that describes this, see Eq. (|38|l and compare with 
the SPA formula in Eq. I)25[l. The approximation is exact 
in both the single particle limit and in the weak tunnel- 
ing limit. Furthermore, we also support our conclusion 
by a perturbation theory in the electron-phonon coupling 
derived in Appendix IbI 



II. MODEL HAMILTONIAN AND 
CURRENT FORMULA 

We consider a model of a single electron level coupled 
to two leads. The single level is coupled to a vibrational 
mode of the molecule through the onsite energy. For 
simplicity we ignore the spin degree of freedom, which 
is not relevant unless Kondo- type effects are important. 
The model Hamiltonian is 



H — H k + H D + H DB + H B + H T , 



(1) 



2 



where 



and 



H k= £fc,« c L c fca, H D =£ dU, 



(2) 



k, a=L,R 



n 2 1 

H B = + 2 m oulxh H DB = Xx d f d, (3) 



(4) 



fc, a=L,R 



Here cc is the oscillator degree of freedom and £o is the 
bare onsite energy. The lead electron energies are given 
by £,ka — £ka — Ha, where /i Q is the chemical potential 
of the lead a. The tunneling amplitude between lead a 
and the molecule is tk a - The amplitudes could in prin- 
ciple also depend on the oscillator position. However, 
in the experimental realizations so far this is a negligi- 
ble dependence and it is omitted here. Because of this, 
we can apply the current formula derived by Wingreen 
and Meir— who expressed the current in terms of the 
non-equilibrium Green's function of the molecule. The 
current through the system is 

1 = I I ^T^TT- K(0 ~ MO] MO, (5) 
n> J 27r 1 L + 1 R 

where the left and right distribution functions are 

n a {e) = n F {e - eV a ), (6) 

and where n F (e) = (e@ e + is the usual Fermi-Dirac 
function. The spectral function, A, is given by 

A(0 = -21mG R (0, (7) 

where the retarded Green's, G R , functions are 

G R (t,t') = -i6(t-t')({d(t),d?(t')}}. (8) 

The widths of the level due to coupling to the leads are 

r a (0 = 27r£ (9) 

k 

The Green's function G R must be calculated in non- 
equilibrium and in the presence of the leads. 

It is often useful to eliminate the linear coupling terms 
of the Hamiltonian Eq. P by a unitary transformation 
similar to one used in the independent boson models 
The transformation is 



H = SHSl S = e - ip ° uU , £ 



X 



(10) 



and with this choice the Hamiltonian transforms into 



H — H k + H B + H D + H T , 



where 



H 



D c 



Sd, 



(11) 



(12) 



Ht = E t ka (c\J^d + cfte-^cj) ■ (13) 



ka 



The retarded Green's function becomes after this trans- 
formation 

G R {t - t') = -i9(t - t)({e i *>® t d(t), d^(t')e- ip o^ e }) s , 

(14) 

where the average of course should be taken with respect 
to H, which is indicated by (•■■)§■ 



III. CALCULATION OF THE 
SPECTRAL FUNCTION 

In this section, we calculate the spectral function, 
A(£) that enters into the current formula, Eq. (J5J. As 
mentioned in introduction, the tunneling broadening of 
phonon assisted side bands has been considered before 
in the case when the presence of the Fermi sea is ig- 
nored(SPA). In this case the model can be solved exactly, 
at least in the so-called wide band limit (WBL)iSiiS For 
latter reference we start by quoting the SPA result 



G 



i?,SPA 



(t)=exp(-tr/2)G£(t), 



(15) 



where G R is the Green's function in the absence of tun- 
neling. We see from this expression that the SPA im- 
plies that all conductance steps are smeared by the same 
amount. 



A. Dyson equation for G 

In the following, we develop a method to calculate the 
broadened Green's function using a truncated equation 
of motion technique. Our starting point is the retarded 
Green's function in Eq. I|14|). We expand it in terms of 
eigenstates of the boson systems 

G R {t) = E Gnn' ,mm' tt) fnn' fmm' i (16) 

G R n ,, mm ,{t) = -i9(t)([(\n)(n'\d)(t),d^\m'){m\}) s , (17) 

where we have defined the overlap function between the 
oscillator states 



/„„, = (n|e*>V>- 



(18) 



An expression for f nn , is given in Eq. (|A12|) . 

We now generate a series of equations of motion, start- 
ing with one for G R n , mm , 

(id t + E nn , - £o)G R n , mm ,(t) = L R n ,(t) 

+ S(t)((dJ\n){m\) s S n , m , + (dU\m')(n'\) s 5 nm ), (19) 
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where E nn , = (n — n')uJ , and L R n , is the contribution 
stemming from the tunneling Hamiltonian. The function 
L R is (see Appendix El for details of the derivations) 

L nn> (*) = E t kafln' ( G k (*)) 



ka.l 



E t kafnfi 



R , (t) 

ridka.,ln . ,mm V / J 



(20) 



where we encounter two new Green's functions 



G f S,, nJ , w (t)=-ifl(t)([(|nXn , |c te )(t),dt| m 'Xm|]> s , 

(21) 

(22) 

The equation of motion for the first one is 

+ 5(i)((c fea dt| nXm | )g ^ m/ + (dt Cfea | m ^|) s J nm ), (23) 

where again L R a n[ is the contribution from the tunnel- 
ing term. Again the tunneling term generates new higher 
order Green's functions, but at this point we truncate it 
using the following physical principles: we neglect cor- 
relations involving lead electrons, which means that we 
decouple terms lead electron operators using the Hartree- 
Fock approximation. Furthermore, we set {c ka <ft) ~ in 
the equation of motion for l|21l) and (|22|l . With these 
approximations (see Appendix lAl for more details) 

Lka,nl{t) ~ ^fc ^ ^ ( n ka) fjn G jl.mm 
3 

+ h a Y,^~^ a ))fi 1 G^. mm (24) 



It is now worth to note that in the case with only one elec- 
tron, which is the SPA, the term (n k ) is exactly zero and 
Eq. I|24|) becomes exact for this case. In the same limit, 
Gn d ka,ni,mm< = 0) aim tne equations are easily solved. 
After setting the result back into Eq. I|16l) we get 

= £ |/ nm ,|» {1 - F n)Nn+ l Nn ' /r (25) 

nn' 

which is nothing but the SPA result in Eq. I|15|) writ- 
ten in the oscillator eigenstate basis. In doing this we 
have furthermore evaluated the last term in Eq. (|19H as 

(<fdt|nXm|)Aw + (d^d\m')(n'\) s S nm = [(1 - n)N n + 
nN n i]5 n , m ,6 nm , where N n = (\n)(n\) is the occupation of 
the n'th oscillator level and n = (cftd) is the average level 
occupation. 

Here we want to go beyond the SPA and we continue 
by looking at the equation of motion for the last Green's 
function in Eq. l(2*2l . We have 



(id t + E nl — £ku)G R d ka,nl,mm>(t) ~ 

Ln d k a ,rd(t)+S(t)(c ka dHn)(m\}6 lm ,. 



The function L R dkanl (t) becomes after doing the same 
line of approximations as was done for L ka nl (t) (see Ap- 
pendix EJ 

" (27) 



L nd ka,nl{t) ~ tka( n ka) ^ , fjn G jl,mm' 



The set of equations now close. After integrating out the 
Green's function involving the lead electrons, we arrive 
at linear equations for the molecule Green's function 

(to + E nn , - eo)G R n , mm ,(uj) = 
{dS\n)(m\) s 6 n , m , + (dU\m){n'\) s 5 nm + 

{■^tn'fnlfjlGfn^mm'^) + ^nlfL' flj G nj,mm' ( W )) ) 

(28) 



E 



where 



n a (0 ) Lu + E ni -£ + irj' 
(29) 

The functions A e ' h are in the case of energy independent 
r a 's (which we assume henceforth) given by 



a 



-^n a (co + E nl ) + ^(co + E nl ) 
Z^n a (-u;-E nl )-^(uj + E nl ) 



where 



(30) 



(31) 



(26) 



where the last term has been added as a regularization 
of the integral at large energies. This is allowed because 
it cancels when adding the two 0j-terms in (|28f) (which is 
easily seen by noting that J2i fnifji = $n,l)- Note that 
■0(e) is a function of e/kT only and that -0(e) = e). 
At large values of e/kT the ^-function has the asymptotic 
form ip(e) ~ — \n(e/kT), while for small e it goes as 
0(e) w ~{e/kT) 2 x 0.213139.... 

The only remaining question is how to evaluate the last 
term in Eq. I|19l) . In absence of tunneling or in the sin- 
gle particle approximation, we can decouple the electron 
and phonon degrees of freedom, such that (dd^ \n)(m\) ss 
(dd' f )N n S nm , where N n = (\n)(n\) is the occupation of 
the n'th oscillator level and n = (d^d) is the average 
level occupation. To leading order in T this is a reason- 
able approximation. However, in the more general case 
we should in principle solve the terms (dd) \n)(m\) and 
(d^ d\m')(n'\) on first right hand of (|28|l self-consistcntly. 
In the general non-equilibrium case, it is, however, not 
possible to calculate these expectation values from the 
retarded Green's function alone. One must therefore in- 
voke new approximations for this purpose. Instead of 
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pursuing that, we concentrate on a case where this is not 
necessary, namely the situation of a strongly asymmetric 
device. This is done in the following section. 

B. Asymmetric transistor, Tl 3> Tr 

In the strongly asymmetric thermal case, things sim- 
plify considerably, because the dot states are in thermal 
equilibrium with one of the leads. For the case 3> Tr, 
the current is 

i= iS^ Tr m ° - n * m al ®> (32) 

where Al is the equilibrium spectral function in equilib- 
rium with the left lead. Therefore, we can use the stan- 
dard equilibrium expression for the self-consistent equa- 
tions 

/oo 
7-\n'.mm'H%H, (33) 
-oo 

/OO 
7-\n',mm'(w)%(-w), 
-oo 27T 

(34) 

where A nn ,, mm ,(w) = i(G R n , <mm , (u) - Gj nn , mml {u)).^ 

We have numerically solved these equations and in 
Fig.^we plot the result for the spectral function, 
for £ Q — 0, and different values of the tunnel broadening. 
It is clearly seen how the Fermi surface effects sharpen 
the lines as compared to the SPA. This effect is more 
pronounced for the first peaks, which is also evident from 
Fig. |21 where we show the spectral function for the off- 
resonance condition. 
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FIG. 1: The broadened tunnel density of states Al(£) for the 
resonance condition, e = 0, g — 2 and two different values 
of r. Top panel: T/uj — 0.5 and lower panel: T/u> = 1. 
The thick solid lines are the numerical solution of the full 
expression, Eq. 1281 , while the thin lines are the single particle 
approximation. Also shown is the approximation expression, 
Eq. 1381 with dashed lines. The approximate expression is in 
fact bare distinguishable from the full result. It is clearly seen 
how the first Frank-Condon step remains sharp even for large 
values of tunnel broadening, in accordance with the estimate 
in Eq. 13911 . In all cases we toke kT/u> — 0.01 



C. Approximate solution of the Dyson equations 

In order to gain more physical understanding of the 
narrowing of the lines seen above and to obtain analyt- 
ical results, we now proceed to solve the equations 
approximately. As is shown in the numerical plots, the 
approximate analytic solution, that is derive below, is in 
fact close to the full solution. 

We look for a solution of Eq. (|28|) near one of the res- 
onances, lu w s + nujQ. Because the Green's function 
G R nn , mm ,{Lu) peaks at u = e - E nn t, the two terms in 
the sum in Eq. 1|28[1 are dominated by, Ej n , « E nn , and 
E n j « E nn , , respectively. Hence we can set j = n in the 
first one and j = n' in second one, which then gives the 
following diagonal equations 



(lu + E nn , - Eq — H nnl {w))G R n , mm ,{to) 
{{l-n)N n +nN n ,}, 



(35) 



where the self-energy is 

JwM - E ( A ^MI/J 2 + A^M|/ in ,| 2 ) ■ (36) 



We have, furthermore, approximated 

(dJ\n)(m\) s 5 n , m , + {dU\m>){n'\) s 5 nm 
~ U„' m '{(l - n)N n + nN n >}, 



(37) 



where N n = cxp(— nf3cu ) and n = (<$d). This approxi- 
mation is valid for not too large T and we have verified 
numerically, in the self-consistent calculation performed 
above, that it is indeed reasonable accurate for T/luq < 1. 
With these approximations, we obtain a closed expression 
for the Green's function G R 



{l-n)N n + nN n , 



lu + E„ 



£0 



\f„ 



(38) 



The result in Eq. (|38|l thus generalizes the SPA approx- 
imation, Eq. (|25l) . by including the leading order effect 
of the Fermi sea. The new self-energy in Eq. has a 
simple physical interpretation: the broadening, which is 
caused by tunneling out and tunneling in processes, can 
only occur if the state in the lead is either empty or oc- 
cupied, respectively. For the non-interacting case, where 
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FIG. 2: Same plot as in Fig. but with e /u) — 1.5 and 
r/u) = 0.5. Note that when compared to the eo = curve 
in Fig.^with the same T (top panel) the first step is broader 
because the blocking of the level broadening gradually disap- 
pears further away from the Fermi energy. 

we in Eq. set E nn i = and use that J^i \f n i\ 2 = 1j 
the spectral functions reduces correctly to an Lorentzian 
with width T. 

We can now in detail study the narrowing of the on- 
resonance line seen in Fig.^ At resonance we have n = 
and we then straightforwardly find 

r ™= |/J2r - (39) 

This result shows that the width of the resonance at 
co = is considerably narrower than the bare tunnel 
broadening would suggest. For large values of g |/ 00 | 2 
becomes smaller and the lines narrows in. However, also 
the weight of the line goes down with |/ 00 | 2 as is evident 
from Eq. (jSHJ- 

D. Comparison with perturbation theory in g 

It is interesting to compare the results where the tun- 
neling is treated approximate but the electron-phonon 
coupling exactly, to the results of a perturbation theory, 
where the tunneling is treated exactly but the electron- 
phonon coupling to lowest order. This situation is solved 
in Appendix [B] and Fig. [3] shows the result of the cal- 
culation. Note that both methods give qualitative the 
same result, namely that the zero energy line is narrowed 
and that spectral function is suppressed between the two 
peaks. 

IV. CONCLUSIONS 

We have presented results for transport through quan- 
tum dots with strong electron-phonon coupling and with 
strong tunneling broadening of the phonon sidebands. 
This is not a solvable problem even when the spin de- 
gree of freedom is ignored, and we resorted to an approx- 
imation which incorporates the Fermi sea, but ignores 




-3-2-10 1 2 3 

CM 



FIG. 3: The perturbative result for the spectral function 
(thick line) compared to the full solution of Eq. 1281 (thin 
line) and the approximate solution, Eq. 13811 . for the sym- 
metric point, e — 0. The perturbative results, which is de- 
rived in Appendix [B] treats the tunneling term exactly but 
the electron-phonon coupling only to linear order. The other 
parameters are F/lo = 0.5, g = 0.25 



correlation effects. The approximation is exact in both 
the single particle and weak tunneling cases. However, 
physical conclusions can be drawn from the approximate 
approach, namely that the tunnel broadening is much 
weaker than expected from a model where the Pauli prin- 
ciple is not incorporated. We have, as mentioned, ne- 
glected the spin degree of freedom which is not expected 
to be important at temperatures above the Kondo tem- 
perature. It is, however, an interesting question to ask 
how the Kondo effect is influenced by the coupling be- 
tween the electron occupancy and the vibrational mode 
of the molecule. 
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APPENDIX A: DERIVATION OF THE 
EQUATIONS OF MOTION 

Here we present the commutators and the methods 
used to derive the set of equations of motion in Section 
IIII Al We write the Hamiltonian in Eq. ljTT|) as H = 
Hq + Ht- In the oscillator eigenstate representation Ht 
is 

Ht = £ t ka (fc^l'Mfw + mViAad) • (Al) 

kot^W 
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Then Eqs. (O,® follow from 

[H ,d\nXn'\] = {E nn , -e )d\n)(n'l 

[H T ,d\n){n'\] = J2 hfw (~dUc ka \l')(n'\S nl - dS c ka \n){lKn>) 



ka.,11' 



The equation of motion in Eq. H23j) follows from the commutators 
[H ,c ka \n}(l\} = (E nl -Z ka )c ka \n)(l\, 

[H T ,C ka \n)(l\] = ^ f k'a' ~fjj> (\n){j'\h a ,k' a ' S jl + c k'a' C ka(\j)( l \ S j' n ~ \ n )(f\ S jlY) d 
k'a'jj' 

+f;r{\f)m n -\n)(j\s jtl )c k/a ,c ka d^. 



(A2) 
(A3) 



(A4) 



(A5) 



At this point we truncate the equations of motion approximating the terms with 3 electron operators by the Hartree- 
Fock decomposition for the lead electrons. We have 



C k'a' C ka d ~ c k' a' { c ka d ^) ~ c ka{ c k' a' d ) ■ 



(A6) 
(A7) 



As mentioned in the main text, we will at this level neglect the terms {c ka d}) and therefore only c k , a ,c ka d contributes 
and gives (c\., a ,c ka )d w 6 ka k , a , (c ka c kol )d. The last approximate sign here means up to order T, which is consistent 
with neglecting (c ka d)). When this is inserted back into l]A5fl . we get 



[H T ,c ka \n){l\] « [-/«' (Wlfy + ( c Lc fea )(b'>ai^„ - \n)(f\5 M ) 

n' 

which then leads to Eq. 124|) . Finally, the last equation of motion in Eq. (|23|) follows from 
[H ,dUc ka \n)(l\] = (E nl - H ka )dUc ka \n)(ll 

[H T ,d^dc k \n){l\]= Y W -fjj'\j)i l \ c k' a ' C ka dS j'n- fjj'\ n )^\ C ka c k' a ' d ^jn 
k'a'jj' 

and after using the same arguments as before, we obtain 

[H T ,dtdc ka \n)(l\] « -J2 t kf jn \3)(l\(4^k a ) d , 



(A8) 

(A9) 
(A10) 

(All) 



which leads to Eq. I]27p. Furthermore, when the oscillator states are chosen real the function / , is 

1 r r~ i i n_ ™'i 

2"2" n!n'! V 2 I v ; 



/ nn' 



(A12) 



APPENDIX B: PERTURBATION THEORY IN g 

In this appendix we derive a perturbative result for the spectral function to first order in the coupling constant 
but to all orders in the tunneling matrix element. We do, however, also assume a constant bare tunneling density of 
states and also take the asymmetric case, <C Tj, = T. We want to compute G R in Eq. (JSJ) and for the perturbative 
calculation it is more convenient to use the Hamiltonian JQ) instead of the transformed Hamiltonian. We will start 
from the equations of motion for the Green's function G R . We use here the notation G R = —i8(t)([y(t),d^]). We 
obtain 



)G?(w) = l + AG&(u;). 



The equation of motion for the function G R d is 



w — £o —i/m 
imuiQ uo — £q 



r<R 
r<R 



(x) 




A=o( w )' 



(Bl) 



(B2) 
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and 

( b " u -tT , ) ( i ) - * ( § ) - « ( : ) g «=»m, <■»> 

where we used that since G^ d in (|B1(1 is multiplied by A, the last term of (|B3|) should be calculate without electron 
phonon coupling. Therefore, for temperatures kT <C luq the following replacements were made in that term Xx 2 d = 
\{x 2 ) d — i£(,d, and Xpxd = X(px) d = =^d. The average value (x) appearing in Eq. i|B2|) follows from (x) — 
—u> 2 {x) — (n d )\/m = 0. Expressing the functions G^ k and G^ k in Eq. (|B3|) in terms of G^ d and G^ d and inserting 
this back into Eqs. I|B2(I and (|B1|I . we obtain 



where N p — lu — £o + P^o + The first four terms in (|B4|) are nothing but the SPA result in Eq. I|25|l expanded to 
lowest order in g. Therefore, it is evident that the last term of Eq. constitutes the correction to the SPA. The 

Green's function G^ k A=n can be evaluated using Wick's theorem 

G&, A =oM = (^>o G £a=oM - (Jc k ) G*^ (u), (B5) 

where 

G^ = oH = tu—^—G*^), (dt Cfc > = i | ^n F (c') (G« A=0 (c/) - G£ A=0 (u/)) • (B6) 

The integrant in Eq. l|B4fl is analytic in the upper half plane of the complex plane, and when performing the 
integral by a contour in the lower half-plane only one term in Eq. (|B6|) contributes. The last term in Eq. (|B4f) . which 
was the correction to the SPA result, then reads 



<^M = 5 r° (B7) 



where 



f°° du' . 1 2cj - Co - . , 

~/(uj)=LJa n F [u})— rpT^— r- (B8) 

The SPA is restored when the energy is either far above or far below the Fermi surface, because the model is then 
similar to a single electron or a single hole. This is clearly seen in the limit |^o|,wo 3> max(r,A:T), where we can 
replace np by either one or zero. For the first cases it is clear that 5G is 0. When np = 1, we have also get zero 
because the integrant is analytic in one half plane. 
Finally, at T — 0, we have 

- ^A'-'w-.i <2, ' n " + ln [,w "° )2 + < T/2 "° )2]) 

(iir6(— ljq — w) +ln|(w +w)/w \) 



2tt(w - £ + w - « r /2) 

u - Co + u) 
2tt(lo — £ — w — *r/2) 



(iir0(ujo -oj) + ln|(w - w)/w |) , 



where 



The zero temperature expression for the SPA correction is seen to diverge at |a>| = luq. This divergence is, however, 
cut-off at finite temperatures. If Fig. [3] we show the result of the perturbation theory compared to the approximate 
spectral function derived in the main text. The plot is at the symmetric point e = Co — S^o = 0- One should note 
that when expanding e to linear order it cancels off the term —guJ /N 2 in Eq. (|B4|1 . 
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